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Abstract 

This  article  presents  a  constitutive  study  of  a  two-phase  composite  based  on  the  authors’ 
previous  paper  [Zhu  H,  Rish  III  JW,  Dass  WC.  Constitutive  relation  for  two-phase  particular 
materials.  I:  elastic  binder.  Computers  and  Geotechnics  1 997 ;20(3):303— 23]  with  inclusions/ 
particles  in  the  composite  behaving  elastically,  but  the  binder/matrix  being  visco-elastic  of 
Maxwell  type.  The  basic-cell  technique  and  binder-contact  laws  are  employed  in  deriving  the 
stress-strain  equations  in  the  integral  representation.  It  is  shown  that  the  overall  stress-strain 
response  of  this  two-phase  composite  can  be  described  by  an  equivalent  multi-axial  visco¬ 
elastic  model  of  four-Maxwell-elements  in  parallel.  The  application  of  this  study  is  aimed  at 
providing  a  mechanics/physics  based  justification  for  employing  the  multi-Maxwell-elements 
model  in  characterizing  asphalt  concrete  and  other  related  issues.  ©  2001  Elsevier  Science 
Ltd.  All  rights  reserved. 

Keywords:  Composite;  Matrix/binder;  Inclusion/particles;  Visco-elastic;  Maxwell  element;  Asphalt  con¬ 
crete 


1.  Introduction 

Compared  with  the  rigidness  of  matrixes  or  binders,  particles  or  inclusions  in 
many  two-phase  engineering  composites  are  much  stiffer.  So,  assuming  particles  or 
inclusions  to  be  elastic  is  a  reasonable  simplification  in  conducting  the  constitutive 
analysis  of  the  composite.  On  the  other  hand,  most  matrixes  or  binders  are  viscous 
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in  nature.  Now  the  question  is,  when  blending  the  two  (matrix/binder  and  particles/ 
inclusions)  together,  what  is  the  overall  constitutive  behavior  of  the  composite? 

There  are  many  studies  done  in  response  to  answering  this  question  in  different 
technical  areas.  Schapery  [3]  and  Hashin  [4]  published  quite  a  few  well-known  arti¬ 
cles  in  the  subject  of  visco-elastic  behaviors  of  composite  materials.  Certain  sands 
exhibit  creep  characteristics.  Murayama  [5]  proposed  a  theoretical  model  to  describe 
it.  Kuhn  [6]  selected  his  PhD  theme  on  the  study  of  sand  creep.  Meegoda  and  Chang 
[7]  attempted  to  incorporate  the  visco-elastic  binder  effect  into  DEM  calculations  in 
simulation  of  asphalt  concrete.  Experimental  studies  of  this  subject  have  also  made 
important  progress,  which  include  the  work  of  Lacerda  [8]  and  Murayama  et.  al.  [9], 
etc.  Zhu  and  Nodes  [10]  presented  a  fabric-tensors  based  approach  in  analyzing  the 
angularity  effect  of  aggregate  in  asphalt  concrete. 

In  the  authors’  previous  article  [1]  (referred  hereinafter  as  part  I),  a  micro¬ 
mechanical  analysis  to  a  two-phase  composite  with  elastic  binder/matrix  and  elastic 
particles/inclusions  has  been  carried  out.  This  article  extends  the  work  done  in  part  I 
in  the  context  that  now  the  binder  is  taken  to  be  a  Maxwell  type  of  visco-elastic 
material.  The  objective  of  this  study  is  to  derive  the  stress-strain  equations  that 
govern  the  mechanical  behavior  of  a  composite  with  elastic  inclusions/particles  and 
a  Maxwell  matrix/binder.  Two  tools  are  employed  to  implement  this  study.  The  first 
one  is  the  basic  cell  method,  and  a  description  of  this  method  is  given  in  part  I.  The 
second  one  is  the  binder-contact  laws  with  Maxwell  type  of  binder  material,  and  the 
details  about  the  laws  can  be  found  in  [2],  There  are  other  types  of  visco  elastic 
binder/matrix  materials,  but  the  reason  to  select  Maxwell  type  in  this  study  is  that 
there  are  quite  a  number  of  investigations  [11,12]  that  are  reported  to  employ  gen¬ 
eralized  Maxwell  models  in  simulating  the  mechanical  behavior  of  asphalt  concrete. 

Since  the  most  steps  in  the  derivation  needed  to  complete  the  current  analysis  are 
identical  to  what  done  in  the  elastic  binder  case,  which  can  be  found  in  part  I,  we 
will  list  only  those  here  in  which  the  viscous  effect  is  presented.  At  the  same  time,  we 
also  refer  to  the  definitions  of  symbols  and  notations  given  in  Nomenclature  in  part 
I  to  both  parts  of  this  work. 


2.  Maxwell  binder  effect 

Fig.  1  shows  a  cross-sectional  view  of  a  two-dimensional  circular  particle  packing 
structure.  The  particles  which  are  all  of  equal  radius,  a,  are  periodically  spaced  in  a 
hexagonal  array  in  a  binder  material  with  the  distance  between  the  centers  of  two 
adjacent  particles  being  d.  Two  packing  directions  known  as  the  closest  packing 
direction  (CPD)  and  the  mid-closest  packing  direction  (Mid-CPD)  can  then  be 
defined.  An  x  y  coordinate  system  is  subsequently  introduced  that  the  x-axis  and  y- 
axis  are  parallel  to  CPD  and  Mid-CPD,  respectively. 

Based  upon  this  packing  structure,  we  select  a  basic  cell  (labeled  A  in  Fig.  1  in  part 
I),  and  will  do  the  analysis  upon  the  configuration  of  this  basic  cell  (see  Fig,  2). 

For  the  elastic  binder,  the  force-displacement  or  stress-strain  is  derived  using  the 
constitutive  relation  which  is  time-independent.  When  the  binder  becomes  visco- 
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elastic,  the  constitutive  relation  now  is  time  dependent.  However,  we  attempt  to 
convert  or  express  such  time-dependent  nature  of  the  constitutive  relation  with  an 
integral  representation  so  symbolically,  the  micro-mechanics  derivation  process  will 


m3 

o  o  o  o  o  o  o 
o  o  o  o  o  o 


Fig.  1.  Sketch  of  the  configuration  for  a  two  phase  material. 


a.  Normal  Forces  b.  Shear  Forces 


Fig.  2.  Basic  cells  under  (a)  normal  forces  and  (b)  shear  forces. 
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be  equivalent  to  that  of  the  elastic  binder  case.  The  relationship  of  pressure-volu¬ 
metric-strain  for  the  Maxwell  binder  can  be  described  by: 


P  ,  P  & 
+  ~  -  Sb 
Ob 


(1-1) 


where,  p  is  the  hydrostatic  pressure,  Kb  and  rjb  are  the  binder  bulk  modulus  and 
viscosity  constant,  respectively. 

Using  the  similar  approach  given  in  Appendix  A  in  part  I,  the  averaged  binder 
pressure  p  in  the  cell  and  the  cell  strains  ex  and  ev  obey  the  following  equations: 

p  =  K*(ex  +  ey)~  <ex  +ey  >*  (1-2) 


where,  <  >  *  represents  an  integral  representation  as  defined  below: 

</>*  =  g*  f  e'**(w,)A*/  dti 
Jo 


(1-3) 


and 


rib 


(1-4) 


where,  Vb  is  the  binder  volume  fraction,  E*  and  K*  are  the  equivalent  bulk  and 
Young’s  modulus  [1]. 

Another  aspect  of  viscous  effect  involves  the  binder-contact  laws.  The  tangential 
and  normal  binder-contact  compliances  for  Maxwell  binder  can  be  generally 
expressed  by: 


Ks&s  —  Fs  H-  §s  Jo^  dt 
K„S„  =  Fn  +  g„f0F„dt 


(l-5a,  b) 


where,  8n  and  8S  are  the  normal  and  tangential  relative  approaches,  Fn  and  Fs  are  the 
normal  and  tangential  contact  forces,  Kn  and  Ks  are  the  normal  and  tangential 
coefficient  of  contact  stiffness.  gn  and  gs  are  defined  as: 


CnbGb 

(Cnb  Cnp)0b 


gs  : 


CsbGb 


(Csb  "1"  Csp)0b 


(l-6a,b) 


and  Cnby  Cnp,  Csb  and  Csp  are  so  called  the  elastic  binder  normal  contact  compliance, 
the  elastic  particle  normal  contact  compliance,  the  elastic  binder  shear  contact 
compliance  and  the  elastic  particle  shear  contact  compliance,  respectively.  Their 
definitions  can  be  found  in  [2]. 

If  we  multiply  the  following  function: 
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eJ>"' 


to  the  both  sides  of  Eq.  (l-5a),  it  gives  rise  to,  after  a  few  steps  of  manipulation,  an 
equivalent  expression  of  Eq.  (l-5a): 

Fs  —  Ks$s~  <  >s  (1-7) 


where,  the  definition  for  <  >s  is: 

<f>s=gs\zU'~")Ksfdh 

Jo 


(1-8) 


It  appears  that  the  term  <8S>S  in  Eq.  (1-7)  represents  the  viscous  portion  of  the 
shear  response  for  the  binder-contact  force-deformation  interaction.  Similarly,  an 
equivalency  to  Eq.  (l-5b)  can  be  expressed  by: 


Fn  =  Kn&n-  <  Sn  >” 


where,  the  definition  for  <  >"  is: 


<f>"=gn 


dti 


(1-9) 


(1-10) 


3.  Basic  cell  analysis 

3.7.  Normal  deformation 


Incorporating  the  Maxwell  viscous  effect  as  stated  in  Eqs.  (1-1)  to  (1-9)  into  the 
derivation  of  the  cell-level  stress-strain  relationships  when  the  cell  is  subjected  to 
normal  forces  (see  the  details  of  derivation  in  part  I),  we  obtain  the  following  two 
integral  equations  between  the  cell  stresses  (crY,  ay)  and  the  cell  strains  (ex,  ey ): 


crx  =  Q\\ex  +  Q\2£y 


Vy  =  Ql2£x  +  Q,22£y 


4-v/3 

V3 


n  V3 

<£x> 

n  3V3 


<  Bx  >"  -- 


<  By  >"  + 


V3 

4 

«  ,V3 


4  <  Ey  >  ”4  <  ex  —  By 


<  By 


>S  ~  <  BX  +  Ey  >* 
>S  ”  <  EX  +  Ey  >* 


(2- la,  b) 


where,  Qn,  Q\2  and  Q22  are  the  elements  in  the  on-axis  stiffness  matrix  [1]. 
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3.2.  Shear  deformation 

As  stated  in  part  I,  the  cell’s  elastic  shear  stiffness  consists  of 
where  Q66 '  and  Q66"  are  non-contact  part  and  contact  part  controlled  shear  stiffness, 
respectively  [1]. 

Following  the  steps  analogous  to  what  given  in  Section  4  in  part  I,  and  at  the  same 
time  incorporating  the  binder’s  visco-elastic  effect,  we  can  obtain  the  cell  shear  stress 
(o^)  —  cell  shear  strain  (exv)  governing  equation: 

a/3  3a/3 

®xy  =  Q66^Sxy  xy  >  jq”  ^xy  *>  <  ^Ex y  >  (2-2) 


where 


<  e.xy  >  =  g 


ejdti,  = 

0  L  3  Tfbs 


(2-3) 


Eqs.  (2-1),  (2-2)  and  (2-3)  are  derived  when  the  cell’s  x  and  y  coordinates  coincide 
with  CPD  and  Mid-CPD  (0  =  0)  [1].  To  put  those  equations  in  a  matrix  form,  we 
have: 


L  -O'  J  fl=0 


fill  Sl2 
fil2  022 
0  0 

<  Ex  >" 


0 
0 
066 


_V3 

"5/3 

1 

O' 

Ey 
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3 

0 

_  y  _ 

4 

0=0 
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1 

<  2 exy  >n 


V3 
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J  0=0 


T  -1 
-l  l 

o  o 
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0 

0.6 


~  <  £x  >s 

"110' 

X 

<  Ey  >S 

— 

1  1  0 

<  lEXy  >S 

0=0 

— 1 

o 

o 

o 

J 

"  <  Sx  >* 

o 

o 

o 

1 _ 

~<EX>' 

X 

<  Ey  >* 

— 

0  0  0 

<  Ey  >' 

<  lEXy  >* 

o 

II 

— 1 

o 

o 
_ 1 

_  <  2e Xy  > 

Je=o 


(2-4) 


4.  Off-axis  analysis 

Eqs.  (2-4)  represent  the  governing  constitutive  equations  for  the  behavior  of  the 
on-axis  cell.  For  the  off-axis  case  (see  the  cell  labeled  B  in  Fig.  1),  when  the  cell 
coordinate  system  is  inclined  at  an  arbitrary  angle  9  to  that  of  the  principal  material 
directions,  the  procedure  of  deriving  the  governing  equations  is  almost  identical  to 
that  corresponding  to  the  elastic  binder  case  (see  part  I,  Section  4).  Hence,  we  here 
only  list  the  final  result: 
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(3-1) 


where  Q\\{9),  Qn{6),  0.22(6),  Q\e{6),  226(0)  and  266(0)  are  the  elements  in  the  off- 
axis  stiffness  matrix.  Their  definition  can  be  found  in  [1].  The  definition  for  5|  and 
Si  are: 

5,  =  cos2(30)  -  sin2(30),  52  =  2cos(30)sin(30)  (3-2) 


The  off-axis  stiffness  matrixes  for  any  given  9  are  presented  in  Eq.  (3-1).  For  a 
packing  structure  which  has  a  statistically  distributed  0  represented  by  a  function 
J{9),  the  averaged  stiffness  responses  then  can  be  evaluated  by  integrating /(0)  with 
both  sides  of  Eq.  (3-1)  over  the  range  O<0<7r/3.  When  J{9)  is  constant,  e.g.  a  uni¬ 
form  distribution,  the  over-all  stress  and  strain  relationship  can  be  easily  integrated 
and  expressed  by: 
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(3-3) 


where  (&x,  axy  and  By)  and  (ex,  exy  and  ey)  are  the  volume  averaged  stresses  and 
strains,  respectively.  2n>  2i2>  Q22  and  266  are  the  elements  in  the  volume  averaged 
stiffness  matrix,  and  their  definition  can  be  found  in  [1], 
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It  is  easily  verified  that  for  all  the  matrixes  in  Eq.  (3-3),  the  identity  remain  true: 
half  value  of  the  difference  between  the  first  and  second  diagonal  elements  is  equal  to 
that  of  the  element  in  first  row  and  second  column.  This  indicates  that,  for  this 
specific  packing  distribution/^)  =  constant,  the  shear  and  normal  deformations  are 
uncoupled,  and  the  stress-strain  relation  defined  in  Eq.  (3-3)  is  isotropic.  This  indi¬ 
cation  appears  consistent  with  the  theme  that  the  uniform  particle  distribution 
(totally  random)  yields  a  solid  state  of  being  overall  isotropic. 


5.  Model  analysis 

All  the  stress-strain  equations  for  the  case  of  on-axis  cell,  off-axis  cell  or  averaged 
response  are  derived  so  far  in  form  of  the  integral  representation  (traditionally, 
those  equations  are  given  in  the  differentiate  format).  It  is  always  desirable  to  see 
what  is  the  visco-elastic  model  those  integral  equations  portray.  For  this  purpose,  a 
model  study  is  conducted  here.  Beginning  with  Eq.  (3-3)  and  by  setting  only  the 
variables  ex  and  ax  non-zero,  now,  the  non-trivial  relationship  between  ex  and  ax 
reads: 


_  — — _  a/35  _  „  \/3  _  _  _  . 

=  Q ii£v - g-  <  sx  >n  — —  <  ex  >s  -  <  ex  >*  -0.5  <  ex  > 


(4-la) 


The  definition  for  Q\\  can  be  found  in  [1]  and  be  expressed  as: 
W\  =  5-^k„ +^-ks  +  k*  +  0.5  Q'66 


(4-lb) 


Consolidating  Eqs.  (4-la)  and  (4-lb)  yields  a  matched  integral  representation 
between  ax  and  ex: 


°x(t)  = 

+ 


V35  K„ 


8 

V3  Ks 


+  K * 


+  0.5  g'66 


£x(t)  -  gn  e  Sn{‘  li)ex(ti)dti 
Jo 


e*(0  -  gs  f 

Jo 


Ut)~g' 


-g'U-t  l)cT-| 


ex(t\)dt  i 


(4-lc) 


It  is  interesting  to  notice  that  the  four  integrals  in  the  right-hand  side  of  Eq.  (4-lc) 
are  of  the  exact  representation  that  characterizes  the  Maxwell  model.  Before  starting 
to  analyze  Eq.  (4-lc),  let  us  introduce  a  model  of  four  Maxwell  elements  arranged  in 
a  parallel  connection  as  shown  in  Fig.  3,  in  which  K\ ,  K2 ,  K3,  and  K4  are  the  stiffness 
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constants,  and  r]\,  Oi*  *?3>  and  t]a  are  viscous  constants  respectively.  It  is  easily  seen 
that  the  following  equations  are  true: 

cr(r)  =  <ri  (r)  +  or2(0  +  0*3  (0  +  0*4(0  (4-2) 


_  0*1  (0  .  0*1  (0  =  0*2(0  0*2(0  _  *3(0  0*3(0  =  0*4(0  <^(0 

K\  &2  m  ^3  m  Ka  oa 


where,  o(t)  and  e(t)  are  the  model  stress  and  strain.  oqW,  cr2(0>  a2(t)  and  cr4(r)  are  the 
stresses  and  strains  defined  in  the  1st,  2nd,  3rd  and  4th  Maxwell  elements  in  Fig.  3, 
respectively. 

Following  the  manipulations  given  in  Eqs.  (l-5)-(l-7),  the  relation  between  e(t) 
and  crt(r)  given  in  Eq.  (4-3)  can  be  expressed  in  form  of  an  integral  representation 
with  the  initial  condition  CTI(r  =  0)  =  A'Ie(r  =  0): 

o-i (0  =  K\e(t)  ”  k\8\  f  e8]{t~h)e(ti)dt\,  g\=—  (4-4) 

Jo  m 

Similarly,  the  relations  between  e(r)  and  cr2(0,  £(0  and  03(0,  and  e(t)  and  cr4(0  can 
be  arrived: 


or2(0  =  K2e(t)  -  K2g2il0e‘g2{t‘ti)B(ti)dtu  g2=  — 
u  m 

03(0  =  *3«(0  -  Kig3jt0emgi(t~tl)e(ti)dtu  gi  -  ^  (4-5) 

0*4(0  =  *4^(0  -  ^AgA^84{t't0s(h)dtu  gA=  — 

u  Oa 

Now,  observation  to  Eqs.  (4-2)  to  (4-5)  easily  proves  that  s(t)  and  a(t)  are  related  by: 
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and  the  comparison  between  Eqs.  (4-1)  and  (4-6)  readily  indicates  that  the  visco¬ 
elastic  model  given  in  Eq.  (4-1)  represents  a  generalized  Maxwell  models  with  four 
Maxwell  elements  in  a  parallel  connection  as  described  in  Fig.  3  with  a  simple  sub¬ 
stitution: 


Ki 


V35  K„ 
8 


Ki  =  K*,  K4  =  O.50'66 


(4-7) 


Ki 

m 


K3  * 

gs*  -  =  g  > 

m 


(4-8) 


The  model  analysis  given  from  Eqs.  (4-1)  to  (4-8)  is  done  to  the  pair  of  ex  and  5X. 
A  similar  conclusion  can  also  be  made  to  the  pair  of  (ey,  ay ),  and  (exy,  axy )  governed 
by  the  second  and  third  equations  in  Eq.  (3-3),  respectively.  In  summary,  the  overall 
uni-axial  visco-elastic  response  of  the  two-phase  composite  can  be  graphically 
depicted  as  shown  in  Fig.  3. 


6.  Asphalt  concrete  simulation 

Asphalt  concrete  (AC)  is  a  multiphase  conglomerate  material  consisting  of 
mineral  aggregates  bound  together  by  a  petrochemical  binder  known  as  asphalt,  and 
is  widely  used  in  civil  and  infrastructure  constructions.  In  this  section,  the  objective 
is  to  examine  a  number  of  subjects  that  relate  to  the  modeling  of  asphalt  concrete 
based  on  the  analyses  and  results  given  in  this  study.  The  starting  point  of  the 
asphalt  concrete  modeling  is  Eq.  (4-1).  Asphalt  concrete  with  good  quality  contains 
typically  3  to  7%  of  air  void.  To  incorporate  this  porosity  effect  in  to  consideration, 
the  modification  is  needed  to  (4-lc).  The  last  two  terms  of  integral  in  the  right-hand 
side  of  Eq.  (4-1)  represent  the  non-contact  portion  of  the  deformation  contribution 
to  its  interaction  with  the  stress.  However,  with  the  presence  of  air  void  in  the 
interstices  of  aggregate  conglomeration,  the  load  transfer  mechanism  within  the 
system  of  this  conglomeration  will  mainly  take  its  path  through  aggregate  contacts. 
As  such,  these  two  terms  are  removed  from  the  stress-strain  equation,  and  corre¬ 
spondingly,  the  equivalent  visco-elastic  model  corresponding  becomes  a  visco-elastic 
model  of  two  Maxwell  elements  in  a  parallel  arrangement  as  shown  in  Fig.  4a.  By 
applying  the  Laplace  transform  technique  to  the  model  given  in  Fig.  4a,  it  can  be 
shown  that  the  visco-elastic  model  in  Fig.  4a  has  an  alternative  representation  which 
is  a  four-element  model  or  Burger’s  model  (see  Fig.  4b).  The  relations  between  the 
parameters  in  Fig.  4a  and  b  are  derived  after  quite  a  few  steps  of  manipulations  as 
follows: 


Ki  =  K\  +  Ki,  oi  =  0\  +  r}2 


gc(Kx  +  K2) 

11,1  (gl  -  gc)(gc  -  g2 )  ’ 


K\  K2  g\K2  +  g2K\ 

gl  = - ,  g2= - .  gc  = 

m  m 


(5-1) 


Ki+K2 


(5-2) 
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Fig.  4.  (a)  Two  paralleled  Maxwell  elements;  (b)  four-element  model. 


KlI  =  Ollgc 


(5-3) 


Using  a  four-element  model  in  simulating  asphalt  concrete  has  been  extensively 
seen  in  public  literature  in  the  area  of  asphalt  concrete  modeling.  In  fact,  it  is  the 
most  employed  approach  in  analyzing  the  visco-elastic  response  of  asphalt  concrete 
specimens.  However,  the  basis  to  use  the  four-element  model  has  been  placed  on  the 
empirical  justification  that  this  model  can  fairly  capture  the  overall  behavior  of 
asphalt  concrete  in  many  aspects,  and  it  is  relatively  simple.  The  study  presented  in 
this  article  actually  provides  a  mechanics/physics  based  justification  that  rationalizes 
the  applicability  of  the  four-element  model  in  its  application  in  describing  asphalt 
concrete. 

Creep  tests  are  often  seen  in  asphalt  concrete  studies.  In  his  work  on  creep  simu¬ 
lation  based  on  the  four-element  model,  Mamlouk  [13]  presents  the  compliance 
relation  as: 


Xt)  = 


e(0 

a-Qii(t) 


K,  +  r]I  +  KII 


0 


(5-4a,  b) 


where,  t  is  the  time  variable;  r  is  retardation  time;  J(t)  is  the  creep  compliance;  e(f)  is 
the  creep  strain;  u(t)  is  the  step  function  and  cr0  is  the  constant  stress. 

Compared  with  E\  and  Eu,  rji  and  On  are  much  more  temperature  dependant.  As 
such,  one  important  subject  in  creep  analysis  is  to  quantify  the  dependence  of  creep 
compliance  on  temperature.  One  well  established  way  to  do  the  quantification  is  to 
introduce  a  set  of  reference  values  for  and  on  at  a  given  temperature.  r]i  and  rju  at 
other  temperatures  then  can  be  linearly  related  to  their  values  at  the  reference  tem¬ 
perature,  which  mathematically  can  be  expressed  by: 

0\  =  mo<ZT\  (5-5) 


r  =  To#T2 


(5-6) 


where,  aT1  and  aT2  are  called  time-temperature  shift  factors.  ^I0  and  r0  are  the  values 
at  the  reference  temperature,  which  is  typically  selected  at  21°C.  Based  on  the 
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argument  that  asphalt  concrete  is  a  thermorheological  simple  material,  Mamlouk 
further  assumes  that 

aj  i  =  aj2  =  aj  (5-7) 

This  assumption  is  also  either  explicitly  or  implicitly  granted  in  most  other  studies 
in  this  subject.  Based  on  Eq.  (5-7),  a  reduced  time  variable  then  can  be  defined: 

$  =  7-  (5-8) 


and  Eq.  (5-4)  now  can  be  written  as: 


J(t)  = 


Bit) 


1  £  1  /,  _i\ 


(5-9) 


The  benefit  of  introducing  £  in  the  creep  compliance  is  that  it  establishes  a  time- 
temperature  correspondence  on  creep  compliance  so  J{t)  determined  by  the  data 
obtained  from  the  test  results  done  at  a  given  temperature  can  be  used  to  estimate 
the  creep  compliance  at  a  different  temperature. 

However,  we  will  demonstrate  here  that  the  assumption  given  in  Eq.  (5-7)  is 
unnecessary,  and  it  is  inherently  valid  without  resorting  to  grant  asphalt  concrete  is 
a  simple  thermorheological  materials.  By  observing  Eqs.  (4-7),  (4-8)  and  (5-1),  it  is 
easily  seen  [also  see  the  definitions  given  in  Eq.  (l-6a,b)]: 


K\  K2 

m  =  m+m=  — i — : 

gn  gs 


{Cnb  "b  ^np) 0b^\  {Csb  *b  Csp)0bK 2 


CnhGb 


CstjGb 


(5-10) 


As  such,  we  introduce  the  shift  factor  for  the  asphalt  binder  only:  Ob  —  Obo  ar  and 
it  will  lead  to  Eq.  (5-5): 


m  = 


(Cnb  H-  Cnp)ribo<3T\K\  ( Csb  -b  CSp)r}tjoaT\K2 


CnbGb 


+  • 


CshGb 


—  OwClT 


(5-11) 


In  an  analogy  and  starting  with  Ob^Obo  t  can  be  expressed  by: 

r]b(K\  +  K2)  0m(K\  -b  Ki)QT 


_  On  _ 

T  KU 


\_ 

gc 


K\  CnbGb  KiCsbGnb 


K\CnbGb  KiCsbGnb 


=  r0flT  (5-12) 


Cnb  "b  Cnn  Csb  ~b  Cnp  Cnb  *b  Cnn  Csb  +  C/ 


-np 


Therefore,  introducing  the  time  shift  factor:  0i  =  rlio  Qr  will  results  in:  r}b—0bo 
and  subsequently:  r  =  to  ax- 
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7.  Conclusions 

The  study  aims  to  develop  a  two-dimensional  constitutive  theory  of  two-phase 
composites  with  the  matrix/binder  being  Maxwell  type  of  visco-elasticity.  The 
stress-strain  relations  are  established  by  using  the  basic  cell  method,  and  the  effect 
of  Maxwell  matrix/binder  is  incorporated  in  the  stress-strain  derivation  by  using  the 
binder  contact  laws.  It  is  noticeable  that  the  derivation  process  is  quite  lengthy  and 
many  integration  and  matrix  manipulations  are  involved.  But,  in  the  end  for  the 
case  of  total  randomness  in  inclusion/particle  packing  distribution,  the  final  result 
appears  that  the  composite  is  isotropic  in  its  homogenization  entity.  This  indicates 
that  introducing  Maxwell  viscosity  in  matrix/binder  does  not  change  the  nature  of 
isotropicity  or  anisotropicity  of  the  composite.  Such  an  indication  is  conceivable.  In 
soil/granular  mechanics,  it  has  been  found  that  complete  uniform  distribution  of 
soil/granular  packing  configuration  will  lead  to  the  isotropic  stress-strain  relation 
defined  in  a  soil/granular  pile  [14],  It  has  also  been  observed  that  a  well  mixed 
asphalt  concrete  exhibits  no  favored  material  orientation  directions,  which  grants 
the  treatment  of  asphalt  concrete  as  an  over-all  isotropic  medium  [15]. 

Most  current  treatment  on  a  visco-elastic  matrix  composite  takes  an  elastic  ana¬ 
logy  approach  that  follows  the  so-called  “corresponding  theory”  [3,4].  In  this 
approach,  the  integral  transform  and  inverse  integral  transform  are  needed,  which 
require  lengthy  mathematical  manipulation  that  the  mathematical  complexity  may 
overwhelm  the  revealing  of  physics/mechanics  for  the  problem.  It  appears  that  the 
study  presents  another  approach  in  investigating  visco-elastic  matrix  composites. 

The  derived  stress-strain  equations  given  in  this  study  are  in  the  integral  repre¬ 
sentation,  and  the  benefit  is  obvious  that  the  visco-elastic  effect  is  concisely  expres¬ 
sed  in  mathematics  and  well  interpreted  mechanically.  The  differentiate  format  of 
the  stress  strain  equations  are  also  available,  however,  very  lengthy  that  the  corre¬ 
spondence  to  Eq.  (3-3)  is  a  set  of  8th  order  of  ordinary  differentiate  equations. 

In  addition,  Eq.  (3-3)  specifies  that  there  are  four  components  that  contribute  the 
visco-elasticity  at  the  composite  level,  and  they  can  be  characterized  by:  (1)  the 
contact  shear  viscous  time  factor  (gs);  (2)  the  contact  normal  viscous  time  factor  (gn); 
(3)  the  non-contact  shear  viscous  time  factor  (g7);  and  (4)  the  non-contact  normal 
viscous  time  factor  (g*).  One  important  characteristic  of  these  four  time  factors  is 
that  they  can  be  analytically  expressed  in  terms  of  the  physical  and  geometrical 
properties  of  matrix/binder  and  particles/inclusions.  In  fact,  this  characteristic  is 
useful  because  it  provides  a  close-form  dependence  of  how  an  individual  parameter 
of  the  binder/matrix  or  inclusions/particles  can  alter  the  behavior  at  the  composite 
level  quantitatively  [16]. 

Using  the  visco-elastic  model  of  multi  Maxwell  elements  in  parallel  connection  to 
simulate  the  behavior  of  asphalt  concrete  is  very  popular  and  can  be  found  in  many 
public  literatures  [15,16,17].  For  example,  Bouldin  et  al.  [16]  showed  that  the 
experimental  strain-time  data  of  an  asphalt  concrete  specimen  under  a  Harversine 
wave  load  fitted  very  well  with  the  curve  computed  by  the  back  analysis  based  on  a 
model  of  two  Maxwell  elements  in  parallel  (Fig.  4a).  In  the  SHRP-A-415,  Moni- 
smith  et  al.  [1 5]  includes  a  section  of  employing  three  Maxwell  elements  in  parallel  to 
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study  the  relaxation  nature  of  asphalt  concrete.  A  four-element  model,  which  is 
equivalent  to  that  of  two  Maxwell  elements  in  parallel,  serves  as  the  basis  in  pre¬ 
dicting  the  creep  compliance  in  Mamlouk’s  work  [13],  The  significance  of  this  study 
is  that  this  study  provides  a  foundation  which  justifies  the  rationality  of  using  the 
model  of  multi  Maxwell  elements  to  represent  asphalt  concrete,  which  is  noticed  by 
the  fact  that  the  four  integrals  in  Eq.  (4-lc)  are  of  the  mathematical  representation 
of  the  Maxwell  model. 

It  has  always  been  an  argument  on  how  many  Maxwell  elements  are  needed  in 
simulating  asphalt  concrete.  The  number  for  being  2,  3  and  4  is  frequently  used. 
From  what  this  study  also  shows  is  that  the  determination  of  this  number  really 
depends  on  the  scope  of  the  simulation  on  the  nature  of  asphalt  concrete.  The  sim¬ 
plest  case  is  what  given  in  Fig.  4a  or  b.  More  importantly,  the  stiffness  and  viscosity 
parameters  in  the  stress-strain  relationships  at  the  composite  level  can  be  analyti¬ 
cally  formulated  on  the  basis  of  physical  and  geometric  properties  of  the  con¬ 
stituents  in  the  composite.  Whereby,  a  parametric  study  can  be  pursued  to  see  how 
the  composite  responds  to  the  change  in  value  of  one  single  parameter  of  physical 
and  geometric  properties.  Such  a  parametric  perturbation  can  be  very  useful  as  an 
analytic  tool  [10]. 
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